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ABSTRACT 

Si x Gei_ x heterostructures on Si substrates have been widely studied due to the maturity of Si 
technology. However, work on Si x Gei_ x heterostructures on Ge substrates has not received much 
attention. A Si x Gei_ x layer on a Si substrate is under compressive strain while Si x Gei_ x on Ge is 
under tensile strain; thus the critical points will behave differently. In order to accurately charac- 
terize high Ge content Si x Gei_ x layers the energy shift algorithm, which is used to calculate alloy 
compositions, has been modified. These results have been used along with variable angle spec- 
troscopic ellipsometry (VASE) measurements to characterize Si x Gei_ x /Ge superlattices grown 
on Ge substrates. The results are found to agree closely with high resolution x-ray diffraction 
measurements made on the same samples. 

The modified energy shift algorithm also allows the VASE analysis to be upgraded in order 
to characterize linearly graded layers. In this work VASE has been used to characterize graded 
Si x Gei_ x layers in terms of the total thickness, and the start and end alloy composition. Results 
are presented for a 1 /zm Si x Gei_ x layer linearly graded in the range 0.5 < x < 1.0. 

INTRODUCTION 

Si x Gei_ x /Si epilayers have recently started to be used for device fabrication. The ability to 
vary the Si composition, x> in the Si x Gei_x layer offers the device designer a wide range of options. 

To date, most work has been focused on SixGei_ x /Si heterostructures grown on Si substrates. 

To our knowledge little work has been done on high Ge content Si x Gei_ x structures grown on 
Ge substrates. The ability to deposit quality Si x Gei_ x /Ge heterostructures will further expand 
the material properties and band offsets available to device engineers. As interest in high Ge 
content Si x Gei_ x materials increases, accurate characterization techniques of such materials and 
structures will be needed. 

Recently, variable angle spectroscopic ellipsometry (VASE) has been shown to be a powerful, 
non-destructive technique for the post-deposition characterization of Si x Gei_ x /Si superlattices 
(SLs) and device structures [1,2|. In these studies VASE was used to determine the layer thick- 
nesses, alloy composition, oxide thickness, sample homogeneity, and number of SL periods. The 
alloy composition was determined by interpolating between the Si x Gei_ x dielectric functions of 
Jellison et ai [3]. However, only the Si rich materials, x > 0.5, were considered. In this current 
work, we have turned our attention to high Ge content Si x Gei_ x /Ge SLs grown on Ge substrates. 

To accomplish this, the Si x Gei_ x interpolation scheme has been modified so that it can be used 
for Si x Gei_ x at all compositions, 0.0 < x < 1.0. 

In order to relieve stress in SixGe^x layers, a graded layer is typically grown on the substrate. 

Graded alloy layers can also be found in the base of heterojunction bipolar transistors (HBT) 
structures. In the past, to characterize structures which contain a graded alloy buffer by VASE, 
the low energy portion of the spectra was ignored, thus ensuring that the graded buffer would not 
affect the analysis. To improve the VASE analysis of structures containing graded buffers and to 
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Figure 1: The Si x Gei_ x dielectric functions of Ref. 3. The dielectric functions for x > 0.5 and 
x < 0.5 were determined from thick Si x Gei_ x layers on Si and Ge substrates respectively. The 
spectra for x = 0.2, 0.47, and 0.98 have been omitted for clarity. 


enable the characterization of graded bases, the analysis software has been updated to account 
for graded alloy layers with an arbitrary alloy profile. 

EXPERIMENTAL DETAILS 

All materials used in this study were grown using a Perkin-Eimer (Model 430S) Si MBE system. 
Deposition of Si and Ge took place in an ultrahigh vacuum (base pressure < 1 x 10“ 10 Torr) from 
dual e-beam evaporation sources onto rotating Si or Ge (100) substrates. After deposition, a 
rotating analyzer ellipsometer was used to record the ellipsometric parameters, {tan 'I', cos A}, 
over the spectral range 300-800 nm. All samples were measured at several angles of incidence 
chosen to increase the sensitivity of the ellipsometric angles to the structure parameters [4]. The 
SL structures were also measured by high- resolution x-ray diffraction (HRXRD) which directly 
determines the thickness and average Si content of one period [5,6]. The thickness and composition 
of the individual SL layers can be calculated from the period, average Si content, and knowledge 
of the shutter opening/closing times from the growth log. 

The energy shift algorithm [7] is used to interpolate between the Si x Gei_ x dielectric functions 
which are given at a number of fixed x values [3,8]. This algorithm requires that the location of 
each critical point (CP) be expressed as a function of the Si content, x . Previous studies of high 
Si content Si x Gei_ x /Si heterostructures [1,2,9-11] have used the functional forms of three CPs 
( E 0 (x ), Ei(x), and Ev (x ) ) in the analysis. The energy shift algorithm has been modified for this 
study in order to account for four CPs by adding the change in the E\ Aj CP with composition. 
The E\ + A i CP has little effect in the energy shift algorithm for relaxed, high Si content Si x Gei_ x 
where the splitting is small, Aj ^ 0.029 at x = 1.0 [12]. The E { + Ai CP function has more of 
an effect for low Si content and strained Si x Gei_ x structures where the splitting is much greater. 

The Si x Gei_x data base of Ref. 3, shown in Fig. 1, has been chosen for this work over the 
data base of Ref. 8 due to its more accurately determined compositions, and wider spectral 
range particularly in the infrared as discussed in Ref. 1. The Si x Gei_ x dielectric functions were 
determined by spectroscopic ellipsometry from thick Si x Gei_ x layers ( d ~ 8/im) on Si substrates 
for x = 98%, 85%, 65%, 47% and on Ge substrates for x - 48%, 28%, 20%, 11%. The E u 
E\ + Aj and E 2 CP functions, shown in Fig. 2, were determined through a CP analysis [12] of the 
Si x Gei_ x dielectric functions. The resolution of the published Si x Ge!_ x data is not sufficient to 
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Figure 2: Ei, E\ + Ai, and E 2 critical points energies calculated from the Si^Gei-* dielectric 
functions of Ref. 3. The CPs determined for Sio.47Geo.53 are displayed as solid symbols and were 
not used to determine the best fit lines. 

resolve the E\ and the E\ + Ai CPs for x > 0.5, so only the E\ CP is considered in this region. 
The dielectric functions for x = 47% and 48% (on Si and Ge substrates, respectively) are quite 
different [3] and lead to inconsistencies in the CP analysis. It can be seen from Fig. 2 that the 
CP energies calculated for x = 47% deviate the greatest from the overall best fit line, and have 
therefore been omitted in the analysis. The CP functions determined from the curves in Fig. 2 
are; 

E x = 2.11 T 1.28a; 

Ei + Ax = 2.28 + 1.13ar (1) 

E 2 = 4.31 - 0.261a: + 0.228a; 2 . 

As E 0 (x) is not an important function in the algorithm, it has been approximated by a linear 
interpolation between the fundamental indirect band gaps of silicon and germanium [13]: 

E 0 = 0.68 + 0.442 . (2) 

Strain dependence was neglected through all the algorithm calculations for two reasons: strain 
effects on the CP were calculated [10, 14] and found to be very small, and no reference dielectric 
functions for strained Si x Gei_ x are available for high Ge concentrations. In addition, measurement 
of such functions is very complicated as the coherently strained calibration samples have to be 
thinner than the critical thickness. 

RESULTS AND DISCUSSION 
SizGei-x/Ge superlattice structures 

Two twenty period Si x Gei_ r /Ge SL structures, labeled A and B, were grown on Ge substrates. 
The target parameters for one period are shown in Table I along with the HRXRD results. 

Each sample was measured using VASE at angles of incidence of 69°, 73°, and 77°. The VASE 
data for Sample A are shown in Fig. 3. To analyze the VASE data a model is generated which 
assumes each layer to be homogeneous and isotropic and all interfaces to be abrupt. In addition, 
each period of the SL is assumed to be identical in terms of both the thicknesses and optical 
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Table I: Comparison of target sample structure with that determined by HRXRD and VASE. 
The period is the sum of the Ge and SiGe layer thicknesses, and x(avg) is the average silicon 
content in one period. The parameters of the VASE analysis are the oxide, Ge, and SiGe layer 
thicknesses, and the silicon content, x, of the SiGe layer. 


Sample 

Source 

Period 

(A) 

x(avg) 

(%) 

d(Ge) 

(A) 

d(SiGe) 

(A) 

X 

(%) 

A 

target 

178.0 

6.2 

128 

50 

22 

A 

HRXRD 

201.1 

8.0 

141.8 

59.3 

26.6 

A 

VASE® 

202.8 

8.4 

127.5 ± 3.9 

75.3 ± 3.5 

22.7 ± 1.6 

B 

target 

201.0 

7.8 

142 

59 

26.6 

B 

HRXRD 

202.1 

8.3 

141.8 

60.3 

27.5 

B 

VASE b 

204.5 

8.4 

126.9 ± 4.7 

77.6 ± 4.4 

22.1 ± 2.0 


a <7 = 0.0137 and d (oxide) = 17.9 A ± 0.4 A. 
h <j = 0.0167 and d(oxide) = 25.5 A ± 0.6 A. 




Figure 3: Comparison of measured VASE data for Si x Gei_ x /Ge SL sample A with the ellipso- 
metric spectra determined from the best fit model. The symbols represent the measured data and 
the lines represent the data calculated from the model. Sample B show similar fitting. 

properties. The native oxide has been simulated using the dielectric function for GeC >2 taken 
from the literature [15]. The results of the VASE analysis are shown in Table I; the period and 
average Si content were calculated from the individual layer thicknesses and x. 

The agreement between the HRXRD and the VASE results are excellent when comparing the 
period and x(avg). The individual layer thicknesses and x, however, differ quite a bit. The {A, 
$} spectra generated from the model are shown in Fig. 3 as dashed lines. The fit can be seen to be 
quite good, especially in the $ spectra. The poor fitting in the long wavelength region of the 69° 
and 73° A spectra is due to the inability of a rotating element ellipsometer to accurately determine 
A values near 0° and 180°. Some small oscillations in A can also be seen in the long wavelength 
region which are not reproduced by the model. In this region the optical penetration depth is 
large and the model is sampling the entire structure. These results suggest that the periods are 
not identical in terms of thickness and/or x(avg) as is assumed in the model. This may in some 
part be responsible for the discrepancies in the individual thicknesses and x. Further work is 
necessary to determine the exact nature of these discrepancies which may also be attributed to 
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Figure 4: Ellipsometric data on graded Si x Gei_ x layers. The experimental data are shown as 
symbols and the best fit model is shown as dashed lines. 

non- uniformities within a single period as a result of a) fluctuations in the Si content of a single 
Si x Gei_ x layer or b) inter facial layers. 


Si x Gei_ x graded layer 


Finally, a 1 /un graded Si x Gei_ x layer was grown on a Si substrate. The layer is linearly 
graded from x — 1.0 to 0.5. The sample was measured by VASE at 70°, 75°, and 77°. The 
measured spectra are shown in Fig. 4. To model the graded layer the total thickness is divided 
into n + 1 (n > 0) layers each of thickness di and Si content Xi. The i-th (0 < i < n) layer 


thickness is given by: 



D/2n 

D/n 


i — 0,n 
1 < i < n 


( 3 ) 


where D is the total layer thickness. The alloy content is assumed to vary linearly between xo 
and x n , 

x = (n - t)xp + ix n ^ 

* n 

In theory this model will approximate a linear continuously graded layer exactly as n -> oo 
although for the graded sample used in this work values of n > 20 yield identical results in the 
VASE analysis. 

The VASE data are modeled by dividing the graded layer into n = 30 layers. The oxide 
thickness is found to be 43.0 A ± 0.4 A, and the unbiased estimator, <r, is 0.0202. The total layer 
thickness, D } is 0.97 /im ± 0.01 /zm which is very close to the nominal values of 1.0 /zm. The 
fitting, shown in Fig. 4, is excellent with the model diverging from the experimental data only in 
the low wavelength region of the V spectra. The alloy composition at the substrate, xo, is held 
constant in the analysis at 100%, the alloy composition at the ambient side of the graded layer 
is 58.7% ± 0.3%. This value differs from the nominal value of 50%. While it is possible that 
the nominal value is incorrect, the discrepancy is most likely due to residual strain in the graded 
layer. Because the graded layer is uncapped, the surface region of the layer is expected to be 
coherently strained. Using relaxed reference dielectric functions in the VASE analysis can cause 
the Si content to be overestimated up to 5% in strained Si x Gei_ x layers on Si [11]. 
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CONCLUSIONS 


The energy shift algorithm has been applied to relaxed Si x Get_ x dielectric functions at dis- 
crete x values covering the complete compositional range. This allows the dielectric function for 
SixGei-x to be generated for all x , 0.0 < x < 1.0. This in turn allows VASE to be used for the 
characterization of high Ge content Si x Gei_ x /Ge heterostructures and also for Si x Gei_ x graded 
alloy layers. VASE results for Si x Gei_ x /Ge SL structures show excellent agreement with HRXRD 
results. VASE has also been successfully used to characterize a linearly graded Si x Gei_ z layer. 
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